Methods for computing a real-time step length and speed of a running or walking individual

ABSTRACT

A method for calculating an estimation of the speed {circumflex over (V)}[n′] of an individual walking along a path, over a time window including: computing a step length  [n′] of the individual, calculating the speed {circumflex over (V)}[n′] using the following formula:
 
 {circumflex over (V)} [ n ′]= [ n ′]× [ n ′].

CROSS-REFERENCE TO RELATED APPLICATION

This application claims priority to European Patent Application No. 18180129.1 filed on Jun. 27, 2018, the entire disclosure of which is hereby incorporated herein by reference.

FIELD OF THE INVENTION

The invention relates to methods for computing a real-time estimation of the step length and the speed of a walking or running individual.

BACKGROUND OF THE INVENTION

Gait speed and step length are important parameters to characterize an individual's daily activities. In sport applications for instance, speed may be used to evaluate athletes and to design personalized training sessions, with a view to improve performance and decrease the risk of injury. In clinical applications, speed is used to evaluate the individual's health, thereby helping doctors diagnose, predict and prevent many diseases such as diabetes, overweight or cardiovascular pathologies.

Global Navigation Satellite System (GNSS) is a basic system widely used to measure an individual's gait speed. GNSS is accurate, and many wearable devices have been designed to embed a GNSS transponder, the measurements of which may be used to compute the individual's gait speed even in real-life conditions. However, there are locations where the GNSS signal is weak or might even be lost due to the lack of satellite coverage, such as indoors, tunnels, near high buildings, narrow valleys. In addition, a GNSS transponder is very power consuming, whereby it is preferable to use it sporadically instead of continuously.

Inertial sensors-based systems are yet another way to measure an individual's gait speed in daily life. While most inertial sensors are used attached to lower limb of the individual's body (such as feet, shanks, thighs), few are used attached to upper limbs (such as trunk and waist).

It is also known to use multiple sensors attached to different parts of the individual's body, the measurements of which are combined to enhance performance of gait speed estimation.

Inertial sensors-based systems estimate gait speed through multiplying gait cadence (i.e. the number of steps per time unit) by step length. Those systems usually compute the cadence by sensing certain moments (e.g. hill-strike, toe-off, mid-swing) during the different gait phases.

In usual inertial sensors-based systems, a step length is computed through geometric modelling of the individual's gait (e.g. based on inverse pendulum movement), double integration of acceleration, or an abstract modelling of the gait based on machine learning.

One main advantage of the inertial-based systems is their possible usage in daily life conditions. Another advantage of the inertial systems is their acceptable accuracy. One main drawback, however, is that those systems need sensor calibration at the beginning of each measurement phase, which requires expertise from the user to wear them on the body. Another drawback is that the location of the sensors on the user's body are often regarded as an inconvenient, especially when a plurality of sensors are required. This might explain their commercial failure.

Only recently were methods introduced to estimate the gait speed of an individual based on wrist-worn inertial sensors.

Most known methods use an abstract modelling, which extracts several features from raw data and then build a regression model to estimate speed. Basic features for the wrist-based speed estimation are energy of the signal, cadence, mean crossing rate, and statistical features like mean, standard deviation, mode, and median of acceleration signal, see e.g. B. Fasel et al., “A wrist sensor and algorithm to determine instantaneous walking cadence and speed in daily life walking,” Med. Biol. Eng. Comput., vol. 55, pp. 1773-1785, 2017. Fasel also proposes intensity of the movement as a feature, which is a multiplication of mean and standard deviation of the wrist acceleration. Fasel yet also proposed the slope of the path (based on pressure sensor) as a feature for speed estimation.

However, even the best methods (including Fasel's) only exploit a little piece of useful information from the physics of hand and body movements. Moreover, the existing methods rely upon a general speed model, which results from an effort to standardize the individual to an average of a large population of humans. While efforts were made to somewhat adapt the general speed model with age, sex, height, and weight parameters, the model remains a standard ignorant of the strategy that each individual might exploit to manage his/her gait speed during physical effort.

SUMMARY OF THE INVENTION

It is an object of the present invention to propose a method for accurate and real-time estimation of the step length of a walking or running individual.

Another object of the invention is to propose a personalized method capable of learning the individual's gait and automatically adapt to his/her activity behaviour.

Another object of the invention is to propose a method that might be implemented within a wrist-wearable device, providing easy-to-use measurements for long-term monitoring of physical activities.

Another object of the invention is to propose a low-consumption computing method to increase autonomy of the wrist-wearable device.

It is therefore proposed, in a first aspect a method for computing an estimation of the step length of an individual over a time window, according to claim 1.

It is proposed, in a second aspect, a method for calculating an estimation of the speed of the individual over the time window, according to claim 10.

It is proposed, in a third aspect, a computer program product implemented on a readable memory connected to a calculation unit, said computer program including the instructions for conducting the method for computing the step length and/or the method for computing the speed as disclosed hereinbefore.

It is proposed, in a fourth aspect, a wrist-wearable device including a calculation unit and a readable memory connected to the readable memory, said readable memory being implemented with the computer program mentioned here above.

BRIEF DESCRIPTION OF THE DRAWINGS

The above and further objects and advantages of the invention will become apparent from the detailed description of preferred embodiments, considered in conjunction with the accompanying drawings.

FIG. 1 is a perspective view of an individual having a wrist-wearable device attached to the wrist.

FIG. 2 is a schematic perspective view of a path followed by an individual during a walking (or running) session.

FIG. 3, FIG. 4 and FIG. 5 are schematic side views of an individual running on a path with no slope (FIG. 3), positive slope (FIG. 4) and negative slope (FIG. 4).

FIG. 6 is a wrist-wearable device embedding a system for computing a real-time estimation of the step length and speed of a walking or running individual.

FIG. 7 is a perspective view showing an individual's hand in movement, with the wrist-wearable device of FIG. 6 attached to it.

FIG. 8 is a diagram summarizing different phases and steps of a method for computing a real-time estimation of the step length and speed of a walking or running individual.

FIG. 9 is a diagram showing two pairs of superposed curves: on the left, a pair of superposed speed curves of a running individual, including a GNSS-built reference curve, and a prediction curve computed with the method of the invention; on the right, the same for a walking individual.

DETAILED DESCRIPTION

On FIG. 1 is schematically drawn an individual on foot 1. The individual 1 is supposed to be healthy and have all his/her members, including legs and arms, the arms having wrists 2. The individual 1 defines a vertical sagittal plane P1 (which is the global symmetry plane of the individual 1) and a vertical front plane P2, which is perpendicular to the sagittal plane P1. As depicted on FIG. 1, a wrist-wearable device or “smartwatch” 3 is attached to at least one of the individual's wrists, here to the right wrist 2.

The smartwatch 3 is configured to monitor activity and infer various parameters of the individual 1 while he/she is moving on foot (i.e. either running or walking) on a path 4, the variations of altitude of which are voluntarily exaggerated on FIG. 2.

Indeed, while some portions of the path 4 may be flat (as shown on FIG. 3), other portions may be inclined and have a slope A which, from the individual's point of view, may be positive (FIG. 4) or negative (FIG. 5).

The individual 1 is moving:

-   -   at a gait speed noted V (expressed in m·s⁻¹),     -   at a cadence noted cad (expressed in s⁻¹ or Hz),     -   with a step length noted SL (expressed in m).

The smartwatch 3 is equipped with:

-   -   a calculation unit 5,     -   a readable memory 6 connected to the calculation unit 5,     -   an accelerometer 7, a barometer 8 and a GNSS transponder 9 all         connected to the calculation unit 5.

The calculation unit 5, readable memory 6, accelerometer 7, barometer 8 and GNSS transponder 9 are embedded in a casing 10 attached to the wrist 2 by means of a strap or bracelet 11.

The readable memory 6 is implemented with a computer program including instructions for conducting a method for computing at least an estimation of the step length over a time window n, noted

[an], of the individual 1 along the path 4, and preferably also an estimation of his/her gait speed over time window n, noted {circumflex over (V)}[n].

The computed step length

[n] and computed speed {circumflex over (V)}[n] are linked to one another by the following equation: {circumflex over (V)}[n]=

[n]×

[n],

where

[n] is an estimated cadence of the individual 1 over time window n.

It should be noted that n is an integer where n ∈ [1, . . . N], N being an integer higher than 1. [1, . . . N] represents a time range starting from a first time window n=1 corresponding to the beginning of the computing to a last time window n=N corresponding to the end of the computing. In a preferred embodiment, time windows are 7-second long with 6-second overlap.

The method uses instant acceleration measures provided by the accelerometer 7 and instant atmospheric pressure measures provided by the barometer 8.

The accelerometer 7 is configured to deliver acceleration measures along three perpendicular axes, namely:

-   -   a first acceleration along an X axis perpendicular to the arm of         the individual 1 and parallel to his/her sagittal plane P1,     -   a second acceleration along a Y axis parallel to the arm of the         individual 1 and parallel to his/her sagittal plane P1,     -   a third acceleration along a Z axis perpendicular to the arm of         the individual 1 and parallel to his/her front plane P2.

For each axis X, Y, Z, z acceleration samples are taken per window, z being an integer higher than 1:

-   -   S_(x) ^(i)[n] is the i-th acceleration sample along the X axis,         delivered by the accelerometer (7) along time window n,     -   S_(y) ^(i)[n] is the i-th acceleration sample along the Y axis,         delivered by the accelerometer (7) along time window n,     -   S_(z) ^(i)[n] is the i-th acceleration sample along the Z axis,         delivered by the accelerometer (7) along time window n.

For convenience, in the rest of the description, the following notations will be used: S _(x)=(S _(x) ^(i))_(i∈[1, . . . z]) S _(y)=(S _(y) ^(i))_(i∈[1, . . . z]) S _(z)=(S _(z) ^(i))_(i∈[1, . . . z])

Likewise, the barometer 8 is configured to deliver z atmospheric pressure samples per window:

-   -   P^(i)[n] is the i-th pressure sample delivered by the barometer         (8) within time window n,     -   P[n] is the average pressure within time window n.

For convenience, in the rest of the description, the following notation will be used: P=(P^(i))_(i∈[1, . . . z]).

The proposed method relies on a machine learning approach based upon the hypothesis that the step length

may be computed using a function, denoted g, of a set of M gait features F₁, F₂, . . . , F_(M) (M being an integer):

=g(F ₁ ,F ₂ , . . . ,F _(m)).

Six gait features (M=6) are considered here, five of which are computed from the acceleration measures S_(x), S_(y), S_(z) delivered by the accelerometer 7, and one of which is computed from the atmospheric pressure measures P delivered by the barometer 8.

F₁[n] is an estimated gait cadence of the individual over time window n: F₁[n]=

[n]. In one preferred embodiment, feature F₁[n]=

[n] is estimated from a known algorithm proposed by Fasel (op. cit. pp. 1775-1778). More precisely, feature F₁[n] over time window n is computed by the steps of:

-   -   computing the norms of the acceleration measures (S_(x)[n],         S_(y)[n], S_(z)[n])=(S_(x) ^(i)[n], S_(y) ^(i)[n], S_(z)         ^(i)[n])_(i∈[1, . . . z]) delivered by the accelerometer 7,     -   applying a Fast Fourier Transform to such computed acceleration         norm,     -   applying a Comb filter to such Fast Fourier Transform, resulting         in a cadence likelihood function,     -   attributing the maximum value of the likelihood function in the         range of frequencies below 2 Hz to F₁[n]=         [n].

F₂[n] is an estimated slope A of the path 4. It is assumed that the step length SL decreases uphill (where slope A is positive, FIG. 4) and increases downhill (where slope A is negative, FIG. 5). In one preferred embodiment, feature F₂[n] over time window n is computed from the atmospheric pressure measures P[n] delivered by the barometer 8. More precisely, a line is fitted to −P[n] using a least square method, and the slope of this line is extracted as the slope of feature F₂[n]. In other words, F₂[n]=A≡−P[n] where the symbol ≡ represents correspondence between A and the linear coefficient of the regression line fitting the variations of −P[n]. This is reflected by the following equation:

${{F_{2}\lbrack n\rbrack} = {- \frac{\sum\limits_{i = 1}^{z}{\left( {i - \overset{\_}{\iota}} \right)\left( {{P^{i}\lbrack n\rbrack} - {\overset{\_}{P}\lbrack n\rbrack}} \right)}}{\sum\limits_{i = 1}^{z}\left( {i - \overset{\_}{\iota}} \right)^{2}}}},$

where l is computed as follows:

$\overset{\_}{\iota} = {\frac{1}{z}{\sum\limits_{i = 1}^{z}{i.}}}$

F₃[n] is an estimated energy of the acceleration measures S_(y) along the Y axis, over time window n. It is assumed there is a neural coupling between the arms and feet movements during gait. Indeed, when the individual increases his step length, the range of arm movement is increased, which leads to an increase in the energy of the acceleration measures S_(y) along the Y axis (the one aligned with the wrist). In one preferred embodiment, feature F₃[n] is computed as the standard deviation of acceleration S_(y) over time window n: F ₃[n]=std(S _(y)[n]), that is to say F ₃[n]=std(S _(y) ^(i)[n])_(i∈[1, . . . z]).

F₄[n] is an estimated mean absolute jerk of wrist acceleration over time window n. During gait periods, especially running, regular impacts are witnessed on the wrist 2 from variations of accelerations, the frequency of which is equal to the step frequency. The mean absolute value of jerk (i.e. the differential of acceleration) provides information about the energy of such impacts. In one preferred embodiment, feature F₄[n] is computed as the mean of the differential of the acceleration measures S_(y) along the Y axis: F ₄[n]=mean(|S _(y) ^(i)[n]−S _(y) ^(i−1)[n]|)_(i∈[1, . . . z]).

F₅[n] is an estimated energy of wrist swing over time window n. It is an estimation of energy of the acceleration measures S_(x) along X axis and S_(z) along Z axis. When the step length is increased during walking, the wrist swing increases, thereby increasing its energy. Feature F₅[n] of wrist swing is computed through the following equation, which provides the standard deviation of acceleration norm of S_(x) along the X axis and S_(z) along the Z axis: F ₅[n]=std(√{square root over (S _(x)[n]² +S _(z)[n]²)}), that is to say F ₅[n]=std(√{square root over (S _(x) ^(i)[n]² +S _(z) ^(i)[n]²)})_(i∈[1, . . . z]).

F₆[n] is an estimated mean of acceleration norm over time window n. During walking, the wrist 2 oscillates around the trunk while the body is moving forward. Feature F₆[n], which denotes a degree of the body movement, is computed as follows: F ₆[n]=mean(√{square root over (S _(x)[n]² +S _(y)[n]² +S _(z)[n]²)}), that is to say F ₆[n]=mean(√{square root over (S _(x) ^(i)[n]² +S _(y) ^(i)[n]² +S _(z) ^(i)[n]²)})_(i∈[1, . . . z]).

The fact that the natures of the movement of wrist 2 during walking and running are different has been taken into consideration. For instant, during walking, the wrist 2 follows a pendulum swing movement that is absent in running.

In order to better fit the model to the type of gait, different combinations of features F₁, F₂, . . . , F₆ are used to implement the method in running and walking, respectively.

More precisely, while F₁, F₂, F₃, F₄ are used for computing function g in a running situation, F₁, F₂, F₄, F₅, F₆ are used for computing said function in a walking situation.

The method for computing the step length

includes two procedures, namely:

-   -   a personalization procedure, including:         -   an offline training phase 100 based upon GNSS data and that             is conducted only once at the beginning of the path 4, and         -   an online training phase 200 which is initiated through the             results provided by the offline training phase and which is             repeated along the path 4,     -   an estimation procedure 300.

Personalisation Procedure

Under the personalisation procedure, the GNSS transponder 9 is ON and delivers instant speeds of the smartwatch 3.

The offline training phase 100 is conducted over an initial time window 0 that does not represent the beginning of the running or walking session but a few seconds later. The offline training phase 100 aims at building an initial model to estimate the step length and the speed of the individual, and includes the steps of:

i) Measuring 102 an initial speed V[0] by the GNSS transponder 9 over the initial time window 0.

ii) Computing 104 an initial step cadence

[0] from acceleration measures S_(x)[0], S_(y)[0], S_(z)[0] delivered by the accelerometer 7 over the initial time window 0.

Initial step cadence

[0] is computed using the steps to compute feature F₁[n], wherein n=0 (

[0]=F₁[0]).

iii) Computing 106 an initial step length SL[0] from the initial speed V[0] and the initial step cadence

[0].

The initial step length SL[0] is inferred from the initial speed V[0] and the initial step cadence

[0] according to the following equation:

${{SL}\lbrack 0\rbrack} = {\frac{V\lbrack 0\rbrack}{\lbrack 0\rbrack}.}$

iv) Computing 108 an initial vector of features X[0] from the acceleration measures S_(x)[0], S_(y)[0], S_(z)[0] and the atmospheric pressure P[0] over the initial time window 0.

The initial vector of features X[0] is computed as follows:

For running: X[0]=[1 F ₁[0] F ₂[0] F ₃[0] F ₄[0]F ₂ ²[0]],

For walking: X[0]=[1 F ₁[0] F ₂[0] F ₄[0] F ₅[0]F ₆[0]].

F₁[0], F₂[0], F₃[0], F₄[0], F₅[0] and F₆[0] are computed using the steps to compute F₁[n], F₂[n], F₃[n], F₄[n], F₅[n] and F₆[n], wherein n=0.

v) Computing 110 an initial coefficient vector β[0] from the initial vector of features X[0] and the initial step length S L[0], following a least square method.

The initial coefficient vector β[0] is computed from the initial vector of features X[0] and the initial step length SL[0]: β[0]=(X[0]^(T) X[0])⁻¹ X[0]^(T) SL[0],

where (X[0]^(T)X[0])⁻¹ is called initial dispersion matrix and is also noted D[0] in a simplified notation.

The previous equation may therefore be written as follows: β[0]=D[0]X[0]^(T) SL[0].

β[0], X[0], D[0], SL[0], V[0] and

[0] are referred to as initial parameters.

The online training phase 200 is performed once the offline training phase 100 is over. The online training phase 200 aims at personalizing the initial model built during the offline training phase 100, and comprises the following steps, sporadically repeated for several time windows n:

i) Measuring 202 an instant speed V[n] by the GNSS transponder 9 over time window n.

ii) Computing 204 an instant cadence

[n] from the acceleration measures S_(x)[n], S_(y)[n], S_(z)[n] delivered by the accelerometer 7 over time window n.

The instant step cadence

[n] is computed using the steps to compute feature F₁. As a reminder, F₁[n]=

[n].

iii) Computing 206 an instant step length SL[n] from the instant speed V[n] and the instant step cadence

[n]

${{SL}\lbrack n\rbrack} = \frac{V\lbrack n\rbrack}{\lbrack n\rbrack}$

iv) Computing 208 an instant vector of features X[n] from the acceleration measures S_(x)[n], S_(y)[n], S_(z)[n] delivered by the accelerometer 7 and atmospheric pressure measures P[n] delivered by the barometer 8 over time window n.

The instant vector of feature X[n] is computed as follows:

For running: X[n]=[1 F ₁[n] F ₂[n] F ₃[n] F ₄[n]F ₂[n]²],

For walking: X[n]=[1 F ₁[n] F ₂[n] F ₄[n] F ₅[n]F ₆[n]].

v) Computing 210 an instant coefficient vector β[n] from the preceding coefficient vector β[n-1] following a recursive least square method.

The instant coefficient vector [n] is computed as follows: β[n]=β[n−1]+D[n]X[n](SL[n]−X[n]^(T)β[n−1]),

where D[n] is a dispersion matrix defined (in theory) by D[n]=(X[n]^(T)X[n])⁻¹. However, to simplify calculation and hence limit memory space, the dispersion matrix D[n] is computed according to the following recursive equation: D[n]=D[n-1](I−X[n](I+X[n]^(T) D[n-1]X[n])⁻¹ X[n]^(T) D[n-1]),

where I denotes the identity matrix.

The recursive computing of equation for computing β[n] starts from β[1]: β[1]=β[0]+D[1]X[1](SL[1]−X[1]^(T)β[0]),

where D[1] is computed using equation for computing D[n] in which the initial dispersion matrix D[0]=(X[0]^(T)X[0])⁻¹ is injected: D[1]=D[0](I−X[1](I+X[1]^(T) D[0]X[1])⁻¹ X[1]^(T) D[0]).

Estimation Procedure

Under the estimation procedure 300, the GNSS transponder 9 is OFF. The estimation procedure aims at estimating the speed, cadence and step length of the individual using only the inertial sensors.

The estimation procedure 300 includes the following steps, for an arbitrary time window n′, n′ ∈ [1, . . . , N]:

i) Computing 302 an instant vector of features X[n′] from the acceleration measures S_(x)[n′], S_(y)[n′], S_(z)[n′] delivered by the accelerometer 7 and atmospheric pressure measures P[n′] delivered by the barometer 8 over time window n′.

ii) Computing 304 an instant step length

[n′] from the instant vector of features X[n′] and the most updated coefficient vector β[n], using the following formula:

[n′]=X[n′]β[n].

iii) Computing 306 an instant cadence

[n′] from the acceleration measures S_(x)[n′], S_(y)[n′], S_(z)[n′] delivered by the accelerometer 7 over time window n′.

iv) Computing 308 an instant speed {circumflex over (V)}[n′] from the instant step length

[n′] and the instant cadence

[n′], using the following formula: {circumflex over (V)}[n]=

[n′]×

[n′].

This method provides an acceptable memory consumption cost and is capable of providing real-time computing of the step length

[n′] and the speed {circumflex over (V)}[n′] over a time window n′. The use of the disclosed recursive least square method for computing the coefficient vectors provides good results without the need of all previous computed data.

Indeed, the method was tested on thirty healthy active volunteers who participated to a 90 minutes outdoor walking and running on a path having various terrain conditions including flat, uphill and downhill.

FIG. 9 shows that the computed speed (in bold line) is close to the actual one (measured through a GNSS, in thin line), whichever the gait (running, on the left; walking, on the right).

This method is inherently user-dependent and hence personalized to each individual. It provides an online learning process where the step length and speed are inferred from instant inputs and from previous computations using a recursive least square model.

The tests conducted on volunteers showed that the error gradually decreases along the path. It proves to have low power consumption, high accuracy, practicality (the method is implemented in a wrist-wearable device) and good ergonomics (no need to calibrate sensors). This method is suitable for a daily use on a long-term basis. 

The invention claimed is:
 1. A method for computing an estimation of a step length

[n′] of an individual running along a path, over a time window n′, n′ being an integer higher than or equal to 1, said individual being equipped with a wrist-wearable device comprising: a calculation circuit, a GNSS transponder configured to deliver instant speeds of the wrist-wearable device, a barometer configured to deliver instant atmospheric pressure measures (P), an accelerometer configured to deliver instant acceleration measures along three perpendicular axes, namely: a first acceleration S_(x) along an X axis perpendicular to a wrist and parallel to a sagittal plane (P1) of the individual; a second acceleration S_(y) along a Y axis parallel to the wrist and parallel to the sagittal plane (P1) of the individual; a third acceleration S_(z) along a Z axis perpendicular to the wrist and parallel to a front plane (P2) of the individual; said method comprising the following steps, carried out by the calculation circuit: computing a vector of features X[n] over a time window n, n being an integer higher than or equal to 1, as follows: X[n]=[1F ₁[n]F ₂[n]F ₃[n]F ₄[n]F ₂ ²[n] where: feature F₁[n] is an instant cadence

[n′] of the individual over time window n, computed from norms of the instant acceleration measures S_(x)[n],S_(y)[n],S_(z)[n] delivered by the accelerometer over time window n, feature F₂[n] is an instant slope of the path computed from the instant atmospheric pressure measures (P[n]) delivered by the barometer over time window n, feature F₃[n] is an instant energy of the second acceleration measures S_(y)[n] delivered by the accelerometer over time window n, and feature F₄[n] is an instant mean absolute jerk of the second acceleration measures S_(y)[n] delivered by the accelerometer over time window n, and computing a coefficient vector β[n] over time window n as follows: β[n]=β[n−1]+D[n]X[n](SL[n]−X[n]^(T)β[n−1]) where: D[n] is a dispersion matrix of the instant vector of features X[n], ${{{SL}\lbrack n\rbrack} = \frac{V\lbrack n\rbrack}{\lbrack n\rbrack}},$ V[n] is a speed of the wrist-wearable device, measured by the GNSS transponder over time window n, β[0] is an initial parameter, computing a vector of features X[n′] over time window n′, as follows: X[n′]=[1F ₁[n′]F ₂[n′]F ₃[n′]F ₄[n′]F ₂ ²[n′]], and computing the step length

[n′] as follows:

[n′]=X[n′]β[n].
 2. The method according to claim 1, wherein feature F₁ [n] is computed by the steps of: computing the norms of the acceleration measures (S _(x) ^(i)[n] S _(y) ^(i)[n] S _(z) ^(i)[n])_(i∈[1, . . . z]) delivered by the accelerometer, applying a Fast Fourier Transform to the computed acceleration norms, applying a Comb filter to such Fast Fourier Transform, resulting in a cadence likelihood function, and attributing the maximum value of the likelihood function in the range of frequencies below 2 Hz to F₁[n].
 3. The method according to claim 1, wherein feature F₂ [n] is computed as follows: ${{F_{2}\lbrack n\rbrack} = {- \frac{\sum\limits_{i = 1}^{z}{\left( {i - \overset{\_}{i}} \right)\left( {{P^{i}\lbrack n\rbrack} - {\overset{\_}{P}\lbrack n\rbrack}} \right)}}{\sum\limits_{i = 1}^{z}\left( {i - \overset{\_}{i}} \right)^{2}}}},$ where: P^(i) [n] is the i-th pressure sample delivered by the barometer within time window k, P[n] is the average pressure within time window n, and ī is computed as follows: $\overset{\_}{i} = {\frac{1}{z}{\sum\limits_{i = 1}^{z}i}}$ where z is number of samples within time window n.
 4. The method according to claim 1, wherein feature F₃ [n] is computed as the standard deviation of second acceleration measures S_(y)[n]: F ₃[n]=std(S _(y)[n]).
 5. The method according to claim 1, wherein feature F₄ [n] is computed as follows: F ₄[n]=mean(|S ^(i) _(y)[n]−S ^(i−1) _(y)[n]) where S^(i) _(y), [n] is the i-th acceleration sample along the Y axis, delivered by the accelerometer within time window n.
 6. The method according to claim 1, wherein D [n] is computed as follows: D[n]=D[n−1](1−X[n]+X[n]^(T) D[n−1]X[n])⁻¹ X[n]^(T) D[n−1]), where D[0]=(X[0]^(T) X[0])⁻¹, and X[0] is an initial parameter.
 7. The method for calculating an estimation of the speed {circumflex over (V)}[n′] of an individual running or walking along a path, over a time window n′, comprising: computing the step length

[n′] of the individual according to the method of claim 1, and calculating the speed {circumflex over (V)}[n′] using the following formula: {circumflex over (V)}[n′]=

[n′]×

[n′].
 8. A computer program product implemented on a non-transitory computer readable medium containing computer instructions stored therein for causing a computer processor to perform the method according to claim
 1. 9. A wrist-wearable device comprising a calculation circuit and a non-transitory computer readable medium connected to the calculation circuit, said non-transitory computer medium containing computer instructions store therein according to claim
 8. 10. A method for computing an estimation of the step length

[n′] of an individual walking along a path, over a time window n′, n′ being an integer higher than or equal to ₁, said individual being equipped with a wrist-wearable device comprising: a calculation circuit, a GNSS transponder configured to deliver instant speeds of the device, a barometer configured to deliver instant atmospheric pressure measures (P), an accelerometer configured to deliver instant acceleration measures along three perpendicular axes, namely: a first acceleration S_(X) along an X axis perpendicular to a wrist and parallel to a sagittal plane (P1) of the individual; a second acceleration S_(y) along a Y axis parallel to the wrist and parallel to the sagittal plane (P1) of the individual; a third acceleration S_(Z) along a Z axis perpendicular to the wrist and parallel to a front plane (P2) of the individual; said method comprising the following steps, carried out by the calculation circuit: computing a vector of features X[n] over a time window n, n being an integer higher than or equal to las follows: X[n]=[1F ₁[n]F ₂[n]F ₄[n]F ₅[n]F ₆[n]] where: feature F₁[n] is an instant cadence

[n′] of the individual over time window n, computed from the norms of the instant acceleration measures S_(x) [n], S_(y)[n], S_(z)[n] delivered by the accelerometer over time window n, feature F₂[n] is an instant slope of the path computed from the instant atmospheric pressure measures (P[n]) delivered by the barometer over time window n, feature F₄[n] is an instant mean absolute jerk of the second acceleration measures S_(y)[n] delivered by the accelerometer over time window n, feature F₅[n] is an estimation of energy of the acceleration measures Sx along X axis and Sz along Z axis over time windows n, and feature F₆[n] is an estimated mean of acceleration norm over time window n, computing a coefficient vector β[n] over time window n as follows: β[n]=β[n−1]+D[n]X[n](SL[n]−X[n]^(T)β[n−1]) where: D[n] is [[the]l a dispersion matrix of the instant vector of features X [n], ${{{SL}\lbrack n\rbrack} = \frac{V\lbrack n\rbrack}{\lbrack n\rbrack}},$ V[n] is a speed of the wrist-wearable device, measured by the GNSS transponder over time window n, and β[0] is an initial parameter, computing a vector of features X[n′] over time window n′, as follows: X[n′]=[1F ₁[n′]F ₂[n′]F ₄[n′]F ₅[n′]F ₆[n′], and computing the step length

[n′] as follows:

[n′]=X[n′]β[n].
 11. The method according to claim 10, wherein feature F₅[n] is computed as follows: F ₅[n]=std(√{square root over (S _(x) ²[n]² +S _(z) ²[n]²)}).
 12. The method according to claim 10, wherein the instant mean F₆[n] of acceleration norm is computed as follows: F ₆[n]=mean(√{square root over (S _(x) ²[n]+S _(y) ²[n]+S _(z) ²[n])}). 